A New Mechanism for Gravitational Wave Emission in Core-Collapse Supernovae 
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We present a new theory for the gravitational wave signatures of core-collapse supernovae. Previ- 
ous studies identified axisymmetric rotating core collapse, core bounce, postbounce convection, and 
anisotropic neutrino emission as the primary processes and phases for the radiation of gravitational 
waves. Our results, which are based on axisymmetric, Newtonian radiation-hydrodynamics super- 
nova simulations (Burrows et al. 2006), indicate that the dominant emission process of gravitational 
waves in core-collapse supernovae may be the oscillations of the protoneutron star core. The oscil- 
lations are predominantly of g-mode character, are excited hundreds of milliseconds after bounce, 
and typically last for several hundred milliseconds. Our results suggest that even nonrotating core- 
collapse supernovae should be visible to current LIGO-class detectors throughout the Galaxy, and 
depending on progenitor structure, possibly out to Megaparsec distances. 

PACS numbers: 04.30.Db 04.40.Dg 97.60.Bw 97.60.Jd 97.60.-s 
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Ever since physicists began to think about gravita- 
tional wave detection, core-collapse supernovae have been 
regarded as prime sources. Traditionally, model calcula- 
tions estimating the supernova wave signature have fo- 
cussed on the gravitational-wave emission from rotating 
iron core collapse and core bounce (see, e.g., Q). Re- 
cent results from stellar evolutionary calculations 0, Q 
and neutron star birth spin estimates 4] indicate that 
presupernova stellar iron cores may rotate much more 
slowly than previously assumed and that the asphericity 
during collapse and bounce due to rotation is not gen- 
erally great enough to produce a sizable time-varying, 
wave-emitting mass-quadrupole moment. In fact, grav- 
itational radiation from large-scale postbounce convec- 
tion and anisotropic neutrino emission are likely to ex- 
ceed the jounce signal of such slowly rotating supernova 
cores 



The first-generation LIGO-class detectors are now op- 
erating at design sensitivity and an international network 
of observatories, including LIGO, GEO600, VIRGO and 
TAMA, is on-line. Gravitational waves detected from a 
supernova can provide us with "live" dynamical informa- 
tion from the supernova core, complementing the super- 
nova neutrino pulse as the only other immediate carriers 
of information from deep inside the star. Using signal- 
processing techniques operating on a large set of theo- 
retical templates, it will be possible to extract supernova 
physics from a sufficiently strong signal j^. 

In recent simulations. Burrows et al. ^ have ob- 
served that protoneutron star (PNS) core g-modes are 
excited by turbulence and by accretion downstreams 
through the unstable and highly-deformed stalled super- 
nova shock (undergo ing the Standing- Accretion- Shock 
Instability [SASI] |1 El El E3) at postbounce times 
of many hundreds of milliseconds. The oscillations damp 
by the emission of strong sound waves and do not ebb un- 



til accretion subsides. In this way the core g-modes act as 
transducers for the conversion of accretion gravitational 
energy into acoustic power that is deposited in the super- 
nova mantle and, as proposed by Burrows et al. 0, may 
be sufficient to drive an explosion. Most easily excited 
is the fundamental £=1 core g-mode, but higher-order 
eigenmodes and, through nonlinear effects, harmonics of 
eigenmodes with complicated spatial structures, emerge 
at later times. 

In this letter, we consider the intruiging possibility 
of the emission of strong gravitational waves from the 
quadrupole spatial components of the PNS core oscilla- 
tions. We obtain new estimates for the gravitational- 
wave signature of core-collapse supernovae. With three 
different presupernova stellar models, we have per- 
formed the longest to date 2D Newtonian radiation- 
hydrodynamics supernova simulations. We find that the 
gravitational waves from the quadrupole components of 
the core oscillations dominate the total wave signature 
in duration, maximum strain, and total energy emis- 
sion by one to several orders of magnitude. We have 
also discovered an approximate progenitor dependence: 
more massive iron cores may experience higher frequency, 
higher amplitude oscillations, and, hence, more energetic 
gravitational-wave emission. 

Method and Initial Models. We carry out our axisym- 
metric calculations with the VULCAN/2D code in the 
mu lti-g roup, flux-limited diffusion approximation H 0, 
The computational grid is comprised of 120 an- 
gular zones on a full hemisphere and 160 logarithmically- 
spaced, centrally-refined radial zones from 30 km out to 
5000 km. The innermost 30 km is covered by a Carte- 
sian region that is deformed to smoothly match the polar 
grid at the transition radius 1]. A resolution test with 
50% more angular and radial zones did not reveal ma- 
jor differences in global integral observables such as the 
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FIG. 1: Dimensionless gravitational- wave strain h+ along the 
equator at a distance of 10 kpc. Note that the range of h+ in 
the lower panel is almost 50 times wider than that of the top 
panel. 



gravitational-wave strain. We employ the Shen equation 
of state [l6]. 

We explore three models in this study. Model sllWW 
is the II-M0 (Zero- Age Main Sequence [ZAMS]) presu- 
pernova model of Woosley & Weaver 17] without rota- 
tion. Model s25WW is nonrotating as well, but is the 
25-M0 progenitor from the same study. Model ml5b6 
corresponds to the IS-M© progenitor model of Heger et 
al. which was evolved with a ID prescription for ro- 
tation and magnetic-field-driven angular momentum re- 
distribution. We map this model onto our 2D grid under 
the assumption of constant rotation on cylinders. It has 
a precollapse ratio of rotational kinetic energy to gravita- 
tional potential energy, p = T/|W|, of -1x10-^%. This 
value is one to two orders of magnitude smaller than in 
previous models (e.g., US El), but yields a PNS con- 
sistent with neutron star birth spin estimates 

We extract gravitational waves from the mass motions 
via the quadrupole formula as described in 

QMS In 

addition, we estimate the gravitational- wave emission by 
anisotropic neutrino radiation with the formalism intro- 
duced by lU and concretized in [2oj. 

Results. Figure depicts the quadrupole gravitational 
wave strain h+ as emitted by mass motions scaled to 
a source distance of 10 kiloparsecs (kpc). In the top 
panel, we superpose the waveforms of models sllWW 
and ml5b6. Despite the presence of some rotation in 
the latter and its greater ZAMS mass, the two models 
have very similar precollapse stellar structures 0, 0, • 
This is reflected in the very similar shapes of their wave- 
forms. Even though sllWW is not rotating, a bounce 
burst strain of ~1.3xl0~^^ (@ 10 kpc) is present in our 
numerical model. The first one to two milliseconds of 
this burst are the imprint of the transition in grid ge- 
ometry from the outer polar to the inner Cartesian grids 
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FIG. 2: Frequency-time evolution of model s25WW's 
gravitational- wave energy spectrum (dEQw/df, [2gI[ ) com- 
puted with a 50 ms sampling interval. 



which generates a time-varying quadrupole moment at 
core bounce. It also induces initial perturbations for vor- 
tical motion in the Ledoux-unstable regions behind the 
expanding shock that sets in almost immediately and 
with a perhaps too fast initial growth rate after core 
bounce. The amount of rotational energy in ml5b6's 
core (at bounce, (3< 0.02% and at the end, ^^-0.08%) is 
too small to have a large influence on the core dynamics 
and, thus, on the waveform, except to slightly stabilize 
the aspherical fluid motion at and shortly after bounce. 
In both models, until about ~250 ms after bounce the 
physical waveform is dominated by convective motions 
in the PNS and in the post-shock region. As the SASI 
0, 0, ^3 ) becomes vigorous and leads to global de- 
formation of the standing shock (see, e.g., j^) the wave 
emission from the post-shock flow increases. 

As described in Burrows et al. 9] the fundamental core 
g-mode (^=l,/~330 Hz in sllWW and ml5b6) is ex- 
cited by turbulence and accretion. It grows strong around 
~400 ms after bounce and starts transferring energy to 
the harmonic at 2/ through nonlinear effects. This is re- 
flected in the rise of si 1 WW's gravitational- wave strain 
around that time. h+ reaches a local maximum, then 
quickly decays to about one- third that amplitude, only 
to pick up again after some tens of milliseconds, rising 
to even higher amplitudes (a maximum of ~7xl0~^^ [@ 
10 kpc]), followed by a quasi-exponential decay with a 
~100 ms e-folding time. We attribute the gravitational- 
wave emission in these two 'humps' to the quadrupole 
spatial component of the 2/ harmonic of the i=l core g- 
mode. si 1 WW's gravitational- wave energy spectrum ex- 
hibits prominent emission in a band around ~650 Hz. A 
frequency analysis shows that the harmonic appears first 
at a frequency of ~590 Hz, which increases over 200 ms 
to a maximum of about 680 Hz, and then continuously 
decreases to ~500 Hz at the end of the simulation. In 
this way, the gravitational- wave emitting component ex- 
actly mirrors the behavior of the i=l g-mode which goes 
through the same phases |2^ . This behavior is qual- 
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itatively consistent with the g-mode frequency evolution 
in PNSs presented in who used general relativity 
and obtained more compact PNSs with higher- frequency 
modes. Interestingly, we find that the time of the first 
lull in the wave emission (~500 ms) coincides with the 
point in time at which the core-oscillation mode reaches 
its maximum frequency and its frequency derivative is 
zero li^. Even though model sllWW begins to explode 
around 550 ms, its core oscillation does not subside im- 
mediately 0. 



Table 1: Model Summary 
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^ at 10 kpc ^ see Eq. (5.3) of ^ see Eq. (12) of Wi 

The lower panel of Fig. ^ displays model s25WW's 
waveform. s25WW's precollapse stellar structure is sig- 
nificantly different from those of sllWW and ml5b6. 
Most importantly, its iron core is more massive (1.92M0 
vs. 1.37M0 and 1.47M0, respectively) and more ex- 
tended. These differences in progenitor structure lead 
to different postbounce evolution for model s25WW. Its 
initial shock radius is significantly smaller and the SASI 
becomes vigorous some 100 ms later than in the two 
other models jl^. Figure El shows a frequency-time plot 
of s25WW's gravitational- wave energy spectrum starting 
at 400 ms after bounce. The first burst of gravitational 
waves, starting at about 500 ms and slowly fading af- 
terwards (Fig. nj, is centered about ~800 Hz. In this 
model, the i=l fiuid mode, which dominates the dynam- 
ics at that time, is centered at ^400 Hz and as for mod- 
els sllWW and ml5b6, we identify the wave-emitting 
component as part of the harmonic at 2/ of the former. 
At ~900 ms after bounce, much stronger waves begin to 
be emitted through the excitation of an i=2 core eigen- 
mode (Fig. [21 at that time / ~950 Hz). It reaches a 
maximum strain of ~5xl0~^^ (@ 10 kpc) and lasts for 
at least 200 ms, emitting a total of close to 10~^ Mq c^ 
(c±1.8xl0^^ erg!) in gravitational waves (Fig.|S|). 

For the three models considered here we find the con- 
tributions of anisotropic neutrino emission to be com- 
pletely negligible compared with those of the core oscilla- 
tions. Extracting the neutrino-fiux anisotropics at an ob- 
server radius of 200 km, well outside the neutrinospheres, 
and extrapolating to 10 kpc distance we find maxi- 
mum strains of 5.5x10"^^, 2.6x10"^^, and 1.3x10"^^ 
for s25WW, sllWW, and ml5b6, respectively. These 
numbers are also up to five times smaller than previ- 
ous estimates obtained with codes that perform neu- 
trino transport along rays, which tends to overemphasize 
anisotropics |3.H.E(]|. while our fiux-limited diffusion ap- 
proach tends to smooth them out slightly. Due to their 
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FIG. 3: Integrated energy emitted in gravitational waves. 

low-frequency character, gravitational waves from neu- 
trinos may not be observable by Earth-based detectors 

Summary. We have derived and analyzed the 
gravitational-wave emissions in long-term Newtonian 
2D radiation-hydrodynamics supernova calculations for 
three different progenitors, one of which includes angular 
momentum consistent with the inferred rotation period 
of young pulsars. We find that the gravitational- wave 
emissions of all models are dominated by the quadrupolar 
components of the core^-mode oscillations first discov- 
ered by Burrows et al. If the core oscillations seen in 
simulations do obtain, their gravitational-wave emission 
will exceed the emission from all previously considered 
emission processes, including rotational core bounce, con- 
vection, and anisotropic neutrino emission. Importantly, 
we point out that the excitation of the core g-modes does 
not depend on the particular explosion mechanism pro- 
posed by Burrows et al. 9], but is likely to be a generic 
phenomenon induced by turbulence and accretion down- 
streams in any suitably delayed explosion mechanism. 

The recent study of Miiller et al. 7] suggests a total 
energy emission of 3xl0~^ MqC^ for a rotating model 
(initially about 3 times as fast as our ml5b6). For our 
nonrotating, low-mass model sllWW we find about 5 
times as much energy emission: 1.6x10"^ Mqc^. For 
the slowly rotating, but otherwise very similar model 
ml5b6, we calculate a total gravitational- wave energy 
of 1.4x10-^ MqC^ (Fig. Eland Table 1). For the ini- 
tial rotation rate and angular momentum distribution 
present in ml5b6, we do not find significant qualita- 
tive or quantitative differences caused by such slow ro- 
tation. Model s25WW, due to its more massive iron 
core, higher postbounce accretion rates, and higher pul- 
sation frequencies and amplitudes, emits an amazing 
8.2xl0~^ Mqc^, almost one tenth of a typical super- 
nova explosion energy. Clearly, radiation-reaction ef- 
fects, not included in our present work, might be sig- 
nificant for this model. In Fig. ^ we show the char- 
acteristic gravitational-wave strain spectra of the three 
models, if located at 10 kpc. /ichar is defined by [i^ as 
hchar{f)=D~^ ^/27^-'^Gc-^dEGw / df , where D is the dis- 
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FIG. 4: Characteristic strain spectra contrasted with initial 
and advanced LIGO (optimal) rms noise curves. 



tance to the source. This is a particularly useful measure, 
since it incorporates the amount of energy radiated in a 
spectral interval df around /. In addition, we show the 
optimal rms noise strain (/irms = yfS{f)^ S{f) being 
the spectral strain sensitivity) of both LIGO I and Ad- 
vanced LIGO 26]. All hchar spectra peak strongly at the 
frequencies identified with the quadrupole components of 
the core oscillations (Fig. [21 between 600 and 1000 Hz, 
likely to be higher when general relativity is included), 
corroborating the narrow-band nature of the emission 
process. Given our results, we conclude that, if the core 
oscillations observed in our simulations are generically 
excited in core-collapse supernovae, even nonrotating su- 
pernovae of small to intermediate progenitor mass should 
be observable by LIGO throughout the Milky Way and 
beyond. Massive progenitors could be detectable out to 
~100 times greater distances and their prolonged and ex- 
tremely energetic core-oscillation wave signature might 
be the generic precursor of stellar-mass black-hole for- 
mation. 

We point out that the work presented here is based 
on simulations in 2D Newtonian gravity and only 
quadrupole wave emission has been considered. Gen- 
eral relativity is likely to increase the frequency of the 
PNS eigenmodes, but is unlikely to lead to qualitative 
differences. Fast rotation might lead to the partial sta- 
bilization of the post-shock convection, might affect the 
growth of core oscillations, and will likely lead to non- 
axisymmetric rotational instabilities for (3 > 8% In 
3D, the temporal and spatial mode and SASI structures 
may change. 
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